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ABSTRACT 

We present a two-zone theory for feedback-regulated star formation in galactic discs, 
consistently connecting the galaxy-averaged star formation law with star formation 
proceeding in giant molecular clouds (GMCs). Our focus is on galaxies with gas surface 
density S g > 100 M pc -2 , where the interstellar medium (ISM) can be assumed to 
be fully molecular. This regime includes most star formation in the Universe and our 
basic framework can be extended to other galaxies. In our theory, the galactic disc 
consists of Toomre-mass GMCs embedded in a volume-filling ISM. Radiation pressure 
on dust disperses GMCs and most supernovae explode in the volume-filling medium. 
A galaxy-averaged star formation law is derived by balancing the momentum input 
from supernova feedback with the vertical gravitational weight of the disc gas. This 
star formation law is in good agreement with observations for a CO conversion factor 
depending continuously on S g . We argue that the galaxy-averaged star formation 
efficiency per free fall time, e^ al , is only a weak function of the efficiency with which 
GMCs convert their gas into stars, e^ c . This is possible because the rate limiting 
step for star formation is the rate at which GMCs form: for large efficiency of star 
formation in GMCs, the Toomre Q parameter obtains a value slightly above unity so 
that the GMC formation rate is consistent with the galaxy-averaged star formation 
law. We contrast our results with other theories of turbulence-regulated star formation 
and discuss predictions of our model. Using a compilation of data from the literature, 
we show that the galaxy-averaged star formation efficiency per free fall time is non- 
universal and increases with increasing gas fraction, as predicted by our model. We also 
predict that the fraction of the disc gas mass in bound GMCs decreases for increasing 
values of the GMC star formation efficiency. This is qualitatively consistent with the 
smooth molecular gas distribution inferred in local ultra-luminous infrared galaxies 
and the small mass fraction in giant clumps in high-redshift galaxies. 

Key words: Galaxies: formation, evolution, ISM, starburst, high-redshift - stars: 
formation 



1 INTRODUCTION 

1.1 Observations of star formation in galactic 
discs and giant molecular clouds 

Galaxies turn their gas into stars remarkably slowly. Aver- 



aged over galactic discs, only e 



gal 



0.01 of the gas mass 



is converted into stars per free fall time in ordinary galax- 
ies (Kennicutt 1998; Genzel et al. 2010). A similar result is 
found when the average is limited to spatial scales > 200 
pc within galaxies (e.g., [Kennicutt et al.||2007 Bigiel et al. 
2008 Blanc et al.|2009"| . However, star formation is not this 
slow down to arbitrary scales. 
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In their study of local molecular clouds, |Evans et al.| 
( |2009[ ) found that the clouds lay a factor of 20 above the 
standard Kennicutt-Schmidt (KS) relation between star for- 
mation rate surface density (E*) and gas surface density 
(Eg). More generally, observations of nearby Galactic clouds 
indicate that such an elevated E* — E g relation is commonly 
found abov e a gas surface density threshold S E ~ 125 Mg 



pc (e.g., Lada et al. 2010| |Heiderman et al. 20101. In- 



terestingly, this gas surface density threshold is comparable 
to the characteristic surface density of molecular clouds in 
the Milky Way and other nearby dwarf and spiral galax- 
ies (e.g., Larson 1981 Solomon et al. 1987 Bolatto et al 



20081. It is thus likely that these departures from the stan- 



dard KS relation on small scales correspond to the transi- 
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tion to gravitationally-bound objects. This interpretation, 
rather than one directly tied to E g , is appealing because the 
disc-averaged E* — E g relation for external galaxies with 
Eg ^> 125 Mq pc~ 2 (including local and high-redshift star- 
burst galaxies) does not show the same elevation^ 

Focusing on the most luminous Galactic free-free 



Murray (20111 inferred the star formation rate per 



free fall time 



GMC _ tsM t 



Mgmc 



(1) 



in giant molecular clouds (GMCs), where tg is the free fall 
time of the GMC, M* is the star formation rate, and Mgmc 
is the total G MC mass (inclu ding stars) . Weighted by ioniz- 
ing luminosity, [Murray | |2oTl| finds (e£ MC ) io n = 0.14-0.240 
This is much larger than the Galactic average, (eff)MW = 
0.006. The difference is significant because most star forma- 
tion in the Milky Way is observed to occur in a relatively 
small number of massive GMCs ( |Williams fc McKee|p97| 



|Murray fc Rahman||20l6"l ). In general, star formation corre 
lates strongly with molecular gas but not atomic gas (e.g., 
Wong fc Blitz|2002|[s"chruba et al.|2011[ > and the mass spec- 



trum of molecular clouds in nearby galaxies is sufficiently 
shallow that the mass is concentrated in the most massive 
objects (e.g., |Williams fc McKee|l997| |Rosolowsky||2005 1. 

Thus, a central problem in understanding what regu- 
lates star formation in galaxies is to reconcile the slowness 
of star formation as observed on galactic scales with the 
higher rate at which (at least some) massive GMCs are in- 
ferred to form stars. Theoretical studies of star formation in 
galaxies have generally not addressed this because they have 
focused either on averages over entire or substantial portions 
of galactic discs (as in the KS law) or on the physics of in- 
dividual GMCs. 



1.2 Theories of galactic star formation 

We review in this section some of the previous results most 
directly relevant to our work. 

|Thompson et ah] ( |2005[ > developed a one-zone model 
of starburst discs supported by radiation pressure on dust. 



Murray et al. ( 2010 1 presented a complementary study of the 



disruption of GMCs by stellar feedback, showing that radia- 
tion pressure is likely the main process by which GMCs are 
disrupted, at least in high-density galaxies where competing 
mechanisms like photoionization and supernovae (SNe) be- 



come inefficient. In particular, Murray et al. (20101 showed 



that in the limit in which GMCs are optically thick to re- 



1 In we compile observations of £g al versus £ g , showing that 
there is no significant trend. 

2 This average does not include a large number of more quiescent 
GMCs, which can have a significantly lower instantaneous e^ MC . 
One interpretation for the wide range of e^ MC values is that it 
is time-dependent, with e^ MC peaking near the end of the life 
of the GMC as it is being disrupted by stellar feedback | |Murray| 
|fc Chang|2012| l. Regardless of the precise interpretation, the ob- 
servations described by |Murray| | |2011[ l demonstrate that GMCs 
hosting ~ 1/3 of the current star formation rate in the Milky Way 
have €g MC well in excess of the galaxy-averaged value. 



processed far infrared (FIR) radiation, a fraction 

RMR M* 



Mgmc 



(2) 



as high as w 0.35 of the initial GMC gas mass is converted 



into stars in its lifetime. Murray et al. ( 2010 1 noted that this 



large integrated GMC efficiency needed to be reconciled with 
galaxy-averaged constraints from the observed KS law. 



Ostriker & Shetty ( 2011[ ) also focused on the starburst 



case. These authors argued that while radiation pressure 
may in fact disrupt GMCs, SNe provide the dominant ver- 
tical support in galactic discs, except in the very innermost 
regions. By balancing the gravitational weight of the disc 
gas with the momentum flux from SN feedback, they de- 
rived a simple expression for the KS law in good agree- 
ment with measurements, provided that the CO intensity 
(^co) to molecular gas surface density (Eh 2 ) conversion fac- 
tor (aco = Eh 2 //co) varies continuously with Jco- Earlier 
work by Ostriker et al. ( 2010 1 considered the complementary 
regime where E g < 100 Mq pc -2 , which we do not consider 
in this work. 



Krumholz et al. ( 2009 1 presented a theory of the KS 
law, covering both E g < 100 M pc" 2 and E g > 100 M Q 
pc -2 , predicting the star formation rate as a function of ob- 
served galaxy properties. This theory is based on the fact 
that star formation proceeds in molecular gas and builds 



on the theory of Krumholz & McKee ( 2005 \ in which the 



slowness of star formation in molecular clouds follows from 
the properties of supersonic turbulence. However, the model 



of Krumholz et al. (20091 leaves some important questions 



open. In particular, it assumes but does not explain how the 
gas is assembled into GMCs with virial parameter a v ir ~ 1. 



Furthermore, the theory of Krumholz & McKee ( 2005 1 does 



not specify how the turbulence is driven and maintained. 



According to the theory of Krumholz et al. (20091, the slow- 



ness of star formation averaged over galaxies derives from its 
slowness within GMCs. The predictions of that theory also 
depend on the assumed fraction of the galaxy mass bound 
into GMCs. In this paper, we derive the GMC mass fraction 
and argue that the star formation rate of galaxies does not 
depend sensitively on how rapidly or efficiently gas turns 
into stars once it is assembled in GMCs. Rather, the critical 
factor is the strength of stellar feedback relative to the grav- 
ity of the galactic disc; the formation rate of GMCs from the 
disc adjusts itself so that the disc-averaged star formation 
law is realized largely independent of sub-GMC physics. 



1.3 This work 

A key element of many star formation theories is the gravita- 
tional instability of galactic discs. Gravitational instability, 
however, can only be a part of a successful theory of star for- 
mation in galaxies. Left unimpeded, gravitational collapse 
would convert all of the gas in a galaxy into stars in about 
one dynamical time, e| al ~ 1, in stark contrast with obser- 
vations. In dense galaxies that are primarily molecular, we 
believe that stellar feedback is the main factor regulating 
star formation. Our goal here is thus to develop a model for 
star formation in GMCs and the resulting galaxy-averaged 
star formation law based on self-regulation by stellar feed- 
back. 

We concentrate on the high gas surface density case 
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(Eg > 100 Mq pc -2 ) to avoid effects related to the conver- 
sion of atomic to molecular gas. Since gas is expected to be 
essentially purely molecular in such discs, we treat the inter- 
stellar medium (ISM) as a single-phase turbulent medium. 



GMCs are disrupted by 
radiation pressure in t GMC 



SNe explode outside GMCs, 
dominate disc turbulence 



This regime includes local merging galaxies (e.g., Downes 
fc Solomon|1998 1, ordinary star- forming galaxies at redshift 



z > 1 (where the star formation rate is elevated at a fixed 
stellar mas s; e.g., |Erb et al.|2006||Daddi et al.|2007||Genzel| 
et al.|2010 1, and high-redshift sub-millimeter galaxies (e.g., 
Tacconi et al.|20"08 |. These systems are among the primary 
targets for new and upcoming observatories sensitive to the 
gas content of galaxies, including the Herschel Space Ob- 
servatory]^] the Jansky Very Large Array (JVLA)Qthe Ata- 
cama Large Millimeter Array (ALMA) jand the Cerro Cha- 
jnantor Atacama Telescope (CCAT)F Since approximately 
half of all stars formed before z ~ 2 and in halos of mass 
comparable to those currently probed by observation (e.g., 
Behroozi et al. 20121, our model in fact applies to most 



star formation in the Universe. The high gas surface den- 
sity regime is also relevant for the inner regions of ordinary 
galaxies like the Milky Way and to the fueling of central 
massive black holes. 

In [J2] we define a two-zone galactic disc model, consist- 
ing of GMCs and a volume- filling inter-cloud medium. In §j3j 
we balance the momentum returned by stellar feedback and 
the weight of the disc to derive an average E* — E g rela- 
tion as a function of the momentum per stellar mass formed 
(P*/m.*) and the Q parameter of the disc. In this section, 
we also derive consistency requirements between the efficien- 
cies of star formation in GMCs and the overall galactic disc, 
e int 1 ° an d e| al - We show that e| al is primarily a function 
of the gas fraction, the circular velocity of the galaxy, and 
the effective P*/m* in the volume-filing medium. There is 
also a dependence on e^ 10 , but it is generally surprisingly 
weak. This is because the disc-averaged Q regulates itself 
to a value that can significantly exceed the threshold for 
marginal stability (Q = 1). This enables the GMC forma- 
tion rate to adjust itself to yield a consistent relationship 
between e| al and e£f c . A comparison with observations in 
S]4] shows that our model agrees well with measurements of 
the galaxy-averaged KS law, the turbulent gas velocity dis- 
persions in galaxies, and the fraction of the disc gas mass 
collapsed in GMCs. We conclude in S|5] Table [T] summarize 
symbols used in this paper. 

While the high surface density limit simplifies the theo- 
retical treatment by avoiding explicit consideration of a mul- 
tiphase atomic and molecular ISM, the main points of this 
paper regarding the interplay between GMCs and the rest of 
the galactic disc likely apply to galaxies with lower surface 
density as well. Similarly, we carry out our numerical calcu- 
lations assuming that supernovae dominate feedback in the 
volume-filling medium and that radiation pressure on dust 
disrupts GMCs, our current best model for galaxies with 
gas surface density 100 < E g < 10 4 Mq pc~ 2 (Appendix 
[A] ). However, many of our results are more general and can 
be extended to other feedback mechanisms (such as pho- 



3 http://herschel.esac.esa.int 

4 http://www.vla.nrao.edu 

5 http://www.almaobservatory.org 
B http://www.ccatobservatory.org 




Radiation pressure dominates 
disc turbulence for S g > 10 4 M Q pc~ 2 

Figure 1. Summary of the two-zone disc models used in this work 
and described in j|2] Star formation proceeds in GMCs, which are 
dispersed by radiation pressure on dust. Most supernovae then 
explode in the volume-filling medium and drive turbulence in it, 
except at £ g > 10 4 Mq pc — 2 , where radiation pressure dominates 
the turbulence in the volume-filling ISM (Appendix [Aj. 



toionization or cosmic rays; e.g. McKee 1989, Socrates et al. 
2008) by appropriate choices for P*/m* and e^t 10 . 



2 TWO-ZONE DISC MODEL 

We consider a two-zone model, consisting of a volume-filling 
background disc in which gravitationally-bound GMCs are 
embedded]^] According to our definition, the self-gravity of 
GMCs exceeds their internal pressure. For any effective 
equation of state in which the pressure increases with av- 
erage density, this implies that the gravity of GMCs also 
exceeds the external pressure acting on them. The GMCs 
are thus treated as collapsed entities that are hydrodynam- 
ically decoupled from the background disc before they are 
dispersed by stellar feedback. Figure [l] shows a schematic of 
our two-zone disc model. 

When comparing our results to observations, it is im- 
portant to distinguish between GMCs as gravitationally- 
bound clouds and the stage in their evolution during which 
they may exist as more compact quasi-virialized objects. 
The mean density of the quasi-virialized state, which fol- 
lows an initial period of gravitational collapse, exceeds the 



7 In the Milky Way, gravitationally-bound clouds roughly coin- 
cide with molecular clouds. In denser galaxies which are the focus 
of this paper, the gas is assumed to be completely molecular. We 
use the term GMC to refer to gravitationally-bound clouds even 
though their chemical composition does not necessarily differ from 
the volume-filling medium. 
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Table 1. Summary of symbols used in this work 



Symbol 



Definition 



Eq. 



e GMC 

,Smc 

e int 
,gal 
% 



K 

h 

M g 
h 

P, »H 
C T 

PT 

A4, M(fc) 
Q, Q(h) 

is, 4 isc 

T 

*GMCi *GMC 

/gmc 

/coll 
a, 

a CO> Xqo 



Star formation efficiency per free fall time of GMCs 1 

Integrated star formation efficiency of GMCs 2 

Star formation efficiency per free fall time in the galactic disc ^] 
Velocity dispersion of isothermal galactic potential, corresponding to circular velocity v c = \/2a 3 

Epicyclic frequency of the disc 5 

Gas disc scale height 5 

Gas mass in the disc 3 

Total and gas mass surface densities 3 

Star formation rate surface density 17 If 8 

Gas mass fraction M g /Mtot H 9 

Mean gas density 3 7 

Turbulent gas velocity dispersion 9~ 

Turbulent gas pressure 1 

Mach number of the turbulence on the outer scale ~ h and for mode k CI 

Toomre stability parameter of the disc and smoothed on scale h B" 

Factor ~ 1 if potential is dominated by stars+dark matter and ~ 1/Q if the gas disc dominates 1 

Free fall time and free fall time evaluated using mean gas density in the disc T2 13 

Momentum returned by stellar feedback in ISM per stellar mass formed T7 TS 

Dimensionless parameters encapsulating ~ 1 uncertainties in — E g relation 19 

GMC lifetime and dimensionless GMC lifetime in units of tf sc 23 [26] 

Fraction of disc gas mass in GMCs 23 

Instantaneous fraction of disc gas mass in gravitationally-unstable density fluctuations 2~S 

Parameters of power-law approximation for "Iqj^q 28 

Conversion from CO intensity to gas surface density and column density 39" 



mean density of the cloud when it first becomes gravitation- 
ally bound. Similarly, the lifetime as a gravitationally-bound 
cloud is longer than that in the quasi-virialized state. In 
this work, we generally define GMC properties in terms of 
their longer existence as gravitationally-bound clouds. This 
frees the model of assumptions regarding how the properties 
of the quasi-virial equilibrium are related to galaxy prop- 
erties. For purely molecular galaxies, gravitational bound- 
edness is also a well-defined criterion for GMCs. However, 
GMCs defined in this way do not in general correspond ex- 
actly to the density peaks traced by molecular gas in the 
Milky Way, which can have an envelope of atomic gas that 
is gravitationally-bound but not fully virialized. 



2.1 Background disc 

Our background disc model is similar to the one introduced 
( [20051 ). 



by Thompson et al. |2005j). The disc is modeled using a 
radius-dependent mean gas density p (including the GMC 
contribution) . As a simple model for galaxies with flat rota- 
tion curve, we assume that the disc is in radial centrifugal 
balance in an isothermal potential with velocity dispersion 
a and angular frequency S7 = y/2a/r. The circular velocity 
is then v c — The total mass enclosed within radius r 

is M tot (r) = 2a 2 r/G and the corresponding surface density 
is 



Etnt — 



TvGr 



(3) 



The gas mass fraction is / g = M s /M to t = E g /E tc .t. 

In general, the vertical component of the hydro- 
static equilibrium equation in cylindrical coordinates is 
(l/p)dp/dz = 9$ g /9z, where $ g is the gravitational po- 
tential. For a thin disc in a spherical gravitational potential, 



this is dp/dz = — pQ 2 z, where p is the effective gas pressure. 
This can be approximated as 

p^4>ph 2 tf, (4) 

where h is the disc scale height and (f> ~ 1 is a constant. In 
the limit in which the self-gravity of the thin disc is domi- 
nant and / g = 1, the solution is similar but with <j> ~ 1/Q- 
Assuming that turbulence dominates the effective gas pres- 
sure^] p w pt ~ per, where ct is the turbulent velocity on 
the scale h. The corresponding Mach number is Mi = ct/cs, 
where c s is the sound speed of the gas. Equation Q implies 
that 



ct 



(5) 



There is strong observational evidence that star-forming 



discs self-regulate to sustain a Toomre ( 1964 1 parameter 



7tGE e 



2<tct 



(6) 



where k = \/4f2 2 + dQ 2 /dlnr = 2a /r is the epicyclic fre- 



quency (e.g., Quirk 1972 Kennicutt 1989 Martin & Ken- 



nicutt 



2001). This can be understood intuitively from the 



8 For a molecular medium at temperature T ~ 100 K, the sound 
speed is only ~ 1 km s , much smaller than obser ved turbu- 
lent velocity dispersions Ct ~ 10 — 100 km s — 1 (e.g., |p~ownes fc| 
|Solomon||1998| |Genzel et al]|2011[ l. In the Milky Way, magnetic 
fields and cosmic rays each contribute comparably to turbulence 
to the total pressure of the ISM (e.g., |Boulares fc 00x11990}. Us- 
ing observations of radio supernova remnants (SNR), |Thompson| 
|et al.| j2009} show that the magnetic pressure is likely a small 
fraction of the total ISM pressure for starbursts with E g > 100 
Mq pc — 2 . Cosmic ray models predict and 7— ray observations in- 
dicate that cosmic rays are also dynamically unimportant for the 
gas surface densities considered here 1 Lacki et al.|20i0 2011 1. 
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fact that galactic discs tend to cool until Q m 1, at which 
point they become gravitationally unstable and form stars. 
Stellar feedback then drives turbulence in the disc, ensuring 
that Q does not drop significantly below unity. We will show 
in this paper that even when stellar feedback can regulate 
the disc to Q ~ 1, significant deviations such that Q > 1 
are expected and are important for regulating the galactic 
star formation rate. Note that in our model, Q is evaluated 
using quantities averaged over the entire disc, including the 
gas in GMCs. 

Real galactic discs consist not only of gas, but also of 
stars and dark matter, each of which has a different surface 
density and velocity dispersion. When the gas turbulent ve- 
locity dispersion is much smaller than the stellar velocity 



dispersion, the results of Rafikov (2001) imply that the ef- 
fective Q criterion for the gas only in equation |6]l is a good 
approximation to the stability of the disc. In the limit in 
which gas and stars have comparable velocity dispersions 
and surface densities, gas and stars contribute similarly to 
the disc instability. This limit may be realized in some cases 
of interest to us, for example in the central regions of galaxy 
mergers ( Downes fc Solomon||1998 1. The stability criterion 
would then be modified by a factor of ~ 2. A more accu- 
rate theory would take this into account, but at the level 
of our analysis we prefer to use the simpler gas criterion in 
equation j6j. 

Using Eg = 2hp and eliminating h using equation (TSp , 
it follows that 



n H 



TrGQr 2 m p 



(7) 



1.7 x 10 4 cm" J Q~ 



200 km s- 



100 pc 



where riH = p/m p . From the definitions of Q and / g , using 
equation we also have the simple relations 

" ° r/ s (8) 



and 



Therefore, 



ct _ Q , 
a 1 



h 
r 



1 ct 
/2(j> a ' 



(9) 



(10) 



2.2 Giant molecular clouds 

As mentioned in the introduction, most star formation in 
galaxies occurs in a relatively small number of massive 
GMCs. We identify the mass of these most massive GMCs 



with the Toomre mass of the disc (e.g., Goldreich & Lynden- 
Bell|1965b, 



(11) 



1.3 x 10 6 M 



^20 pcy vio 3 m q pc- 2 

where we have scaled the parameters to values characteristic 
of local ULIRGs. The Toomre-mass GMCs are initially self- 
gravitating but are eventually dispersed by feedback from 



massive stars forming in them (e.g., Matzner 2002; Murray 
et al. 20101. We assume that all stars form in GMCs of 



Toomre mass and denote the fraction of the disc gas mass 
stored in gravitationally-bound GMCs at any given time by 
/gmc: Mgmc = /GMcAfg. 

In general, the free fall time in a region of mean density 

p is 



tff 



371 



32Gp 



(12) 



-1/2 



Since the mean density of a GMC in quasi-virial equilibrium 
exceeds the mean density in the disc, the free fall time of 
the GMC in that state is shorter than the free fall time in 
the disc. In the Milky Way, quasi-virialized GMCs are over- 
dense with respect to the disc by a factor ~ 30, so that their 
internal free fall time is shorter than the free fall time in the 
disc by a factor ~ 5. The lifetime of quasi-virialized GMCs 
before they are dispersed by feedback is likely a few internal 
free fall times (e.g., [Murray fc Cha ng 2012), comparable to 
or less than the free fall time at mean disc density, tg lsc . Since 



the disc (for Q ~ 1), the GMC lifetime as a gravitationally- 
bound entity £gmc ~ iff sc . 

Using equation Q for the mean gas density in the disc, 



.disc 
Iff = 



3tt 2 Q 



32 x 2V20V2 a 
- Q 1 ' 2 



(13) 



4 x 10° yr 



<fjV 4 ^ 100 pc/ V 200 km s 



The short dynamical times in dense starburst discs imply 
that massive stars can outlive their parent GMC. We ex- 
pand on this point in Appendix [X] in the context of which 
feedback processes disrupt GMCs and which processes drive 
turbulence in the volume-filling ISM. 



3 CONNECTING DISC AND GMC STAR 
FORMATION 

We now connect star formation in GMCs and the disc- 
averaged star formation law. We first derive a general ex- 
pression for the disc-averaged KS law based on vertical hy- 



drostatic equilibrium in the disc ((3.11, discuss the relation- 



ship to the turbulent gas velocity dispersion (i 3.2 1, and show 



how the implied disc-averaged star formation efficiency is re- 
lated to the integrated star formation efficiency and lifetime 



of GMCs (p3l. 

We express our results here in terms of a fiducial value 
for P*/m* = 3,000 km s _1 (the effective momentum in- 
jected by stellar feedback per stellar mass formed) appropri- 
ate for SN feedback under typical conditions. We motivate 
this choice and discuss how P*/m* depends on the ambient 
conditions (including the importance of radiation pressure) 
further in Appendix [A"] 



3.1 The disc-averaged star formation law 

The general requirement that the turbulent pressure bal- 
ances the gravitational weight of the overlying gas in the 
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disc is 



Pt = 



7rG£ 6 
2 3 / 2 



(14) 



This expression follows from hydrostatic equilibrium and the 
definition of Q, and the pre-factor depends on the assump- 
tion of an isothermal potential. Equation ( 14 1 holds for any 
value of Q or / g . A similar equilibrium relation was derived 



by |Thompson et al.| pOOH) and |Qstriker fc Shetty| ( |2011 1. 

In general, the total turbulent pressure results from a 
combination of stellar feedback and all other processes that 
drive turbulence: 



PT = PT,* + PT.othe 



(15) 



The second term, Pt, other, may include contributions from 



the magneto-rotational instability (MRI; e.g., Piontek & Os- 
triker|2005 I, the thermal instability (e.g., Kritsuk & Norman 
2002 1, the release of gravitational potential energy as gas ac- 



cretes onto the galaxy from the intergalactic medium (e.g., 
Genel et al.|2012 1 or as gas is transported inward by internal 
galactic torques (e.g., Krumholz fc Burkert|2010" l. 

We focus on the case in which stellar feedback domi- 
nates turbulence, pt ~ Pt,*- Then pt,* can be derived by 
equating the rate of energy injection per unit volume, ein,*, 
with the total turbulent energy dissipation rate, edi ss - The 
energy injection rate ei n ,* ~ E*(P + /m*)n*/2, where is 
a velocity term used to convert momentum to kinetic en- 
ergy injection. If isolated SNe dominate the feedback, is 
the supernova remnant (SNR) velocity at the stage used to 
evaluate P^/m*. We define this stage such that d 4 = ct, 
which for SNRs corresponds to when the remnants effec- 
tively merge with the ISM. In general, P*/m* can also be 
determined by other processes, such as photoionization in 
lower-E g galaxies (e.g., |McKee|1989|[Qstriker et al.|2010[ ) or 
radiation pressure on dust in higher-£ g galaxies (Thompson 
et al. 2005 and Appendix |A|) . 

Numerical simulations show that turbulence dissipates 
in approximately one flow crossing time, 

L 



ct 



(16) 



where 7 « 1 and L is the size of the largest eddies (e.g., 
Stone et al. 1998 Mac Low 1999 1. This expression is valid for 



hydrodynamic and magnetohydrodynamic turbulence, both 
subsonic and supersonic. The turbulent energy dissipation 



rate is thus edi s 
for pt, 



pc%: 



pc T /idiss = P c t/{"iL) an d we can solve 



„ _ fpQ- - /w)/ft7 yi 

Pt,* ~ 



(17) 



In the above equation, we defined fh = L/h and introduced 
a term /p(l — /w)/4, where the factor of 1/4 accounts for 
cancelation of momentum in the disc plane and fp param- 
eterizes the uncertainty in it (Ostriker & Shetty 2011 1. We 
also defined / w as the fraction of the input momentum that 
is lost to a galactic wind, rather than contributing directly 
to the vertical pressure support in the disc. 

The dependence on fh = L/h in equation (171 illus- 



trates the dependence on the details of how turbulence is 
driven. For example, if turbulence is driven by individual 
SNe and the ambient density is very large, then it is possi- 
ble that L <C h. In this limit the turbulent energy is ef- 
ficiently radiated away and SNe are relatively inefficient. 



On the other hand, SNe could cluster and merge in super- 
bubbles of scale ~ h before achieving pressure equilibrium 
with the ambient ISM. 

Equating equations ( |14[ ) and ( |17[ ), we obtain the star 
formation law 



(18) 



13 M yr _1 kpc~ 



v 3,000 km s- 1 J V 103M eP c_2 / ' 
where 

T=fp(l-U)f hl (19) 
encapsulates uncertain factors of order unity. 



3.2 Relation between turbulent gas velocity 
dispersion and star formation efficiency 

We can also write the star formation law by defining a disc- 
averaged star formation efficiency per free fall time, e| al , as 
in the usual KS law: 



(20) 



Using the definition of tg 1BC (eq. 12 1 and expressing the mean 
density p using the vertical balance equation ( 14 1 , the star 
formation law becomes 



,7/4 



GQ 



73 



ct 



(21) 



We can further use equation ( 18 1 to eliminate E g and solve 
for the turbulent velocity: 



,1/4 



c T = 



(22) 



6.6 km s 



Ql/201/2 I o.oi 



3,000 km s- 



This result implies that the turbulent velocity disper- 
sion is constant at fixed e| al , even if the supernova rate 
varies by orders of magnitude with the star formation rate, 
in agreement with the numerical simulations of Joung et al. 
( 2009 ) . In j |4.4[ we show that the elevated turbulent velocity 
dispersions observed in local ULIRGs and in high-redshift 
star-forming galaxies can be explained by the dependence 
of e| al on the product / g cr, which we discuss next. 



3.3 Consistency with star formation in GMCs 

Consider the total mass of stars formed in a time interval of 
duration 4 isc , M*(4 isc ). This can be written in two ways: 

M*(4 isc ) = (P-^-) M g / G Mce£f C (GMC view) (23) 

V'GMC / 



and 

AL(4 isc ) = e| al M g 

The right-hand side of equation ( |23[ ) is simply the number of 
GMC generations in one tg BC times the stellar mass formed 



(disc average view). (24) 
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Figure 2. Equilibrium solutions for the disc-averaged Toomre Q parameter, fraction of the gas mass in gravitationally-bound GMCs 
(/gmc)i an d the disc-averaged star formation efficiency per free fall time (e^ al ) as a function of the integrated efficiency with which 
GMCs convert their gas into stars Each panel shows the curves for different Mach numbers M of the molecular ISM. The left 

panel is representative of a low-mass or gas-poor galaxy with f g <r = 20 km s _1 and the panel on the right is illustrative of a massive 
gas-rich system (such as an ULIRG or a high-redshift star-forming galaxy) with f g cr = 100 km s _1 . The disc-averaged star formation 
efficiency per free fall time increases slowly with the integrated GMC star formation efficiency, while /gmc decreases. This is realized by 
the global Q parameter increasing above the classical stability threshold of unity, so that GMCs form only where turbulent fluctuations 
cause the self-gravity of the gas to exceed its turbulent and rotational support. The se examples assume P*/m* = 3,000 km s _1 , T = 2 
(shown to yield a good fit to the observed disc-averaged X* — S g relation in [4.2 I, <j> = 1, and Jgmc = 1- The equilibrium solutions 
were obtained using equation | |C7| for the collapsed fraction as a function of Q, with values (a, b, \) = (1.5, 0.5, 1) for the parameters 
defined in Appendix C The curves start at ej^"^ = 0.02 in each panel. 



in GMCs in one generation, while equation ( 24 1 follows di- 
rectly from the definition of e| al [^]As before, /gmc and £gmc 
are the disc gas mass fraction in GMCs and the total GMC 
lifetime, where GMCs are defined as gravitationally-bound 
clouds (see Sj2|. Since we assume all star formation occurs 
in GMCs, we can equate these two expressions to find 



gal _ /gmc gmc 



es = - 



£gmc 



(25) 



9 In equations j23[ l and I I24J , the time interval is arbitrary and 
tg lsc was chosen to simplify the algebra that follows. It is assumed 
that the GMC lifetime is at least one free fall time at the den- 
sity threshold for gravitational instability, which ensures that the 
turbulence re-arranges itself rapidly enough that the number of 
GMC generations in a given time interval is determined by the 
GMC lifetime, rather than the time scale for new density fluctu- 
ations to begin gravitational collapse to a GMC. 



where 



tGMC 



\ ^disc / 
V % / 



(26) 



is the GMC li fetim e as a fraction of the disc free fall time. 
We argued in \ 2.2 that £gmc ~ 1- By definition, 



/gmc 



7kfg/4 isc 



(27) 



£gmc 

The numerator is the rate at which gas is processed 
by GMCs, i.e. either turned into stars or returned to the 
ISM. The denominator is proportional to the rate at which 
gas from the disc is incorporated into bound GMCs, Mgmc- 
More precisely, the GMC formation rate is determined by 
the rate at which turbulent fluctuations become gravitation- 
ally unstable, 

Mg/ co u 



GMC 



Xtflo 



(28) 



where / co n is the fraction of the gas mass in the disc that 
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Figure 3. Equilibrium solutions for the disc-averaged Toomre Q parameter, fraction of the gas mass in gravitationally-bound GMCs 
(/gmc)i arl d the disc-averaged star formation efficiency per free fall time (e|- al ) as a function of the galaxy gas surface density (E g ). Each 
panel shows the curves for different Mach numbers Ai of the molecular ISM. The left panel is representative of a low-mass or gas-poor 
galaxy with f g a = 20 km s _1 and the panel on the right is illustrative of a massive gas-rich system (such as an ULIRG or a high-redshift 
star-forming galaxy) with / K cr = 1 00 km s" 1 . We assume a simple model of GMC disruption by radiation pressure on dust to evaluate 
e?J^ c as a function of E g (eq. |38| . The model predicts that Q and e|. al increase modestly with Eg, while /gmc decreases. The flattening 
at Eg = 1,600 Mq pc~ 2 is caused by the saturation of e^" 10 in the optically thick limit. These examples assume P t /m t = 3,000 km 
S — 1 , T = 2 (shown to yield a good fit to the observed disc-averaged E* — E g relation in [4.2 i, </> = 1, and tGMC = 1- The equilibrium 
solutions were obtained using equation | |C7| for the collapsed fraction as a function of Q, with values (a, b, \) = (1-5, 0.5, 1) for the 
parameters defined in Appendix [C] 



is unstable to gravitational collapse at any time. We have 
identified the time scale for the turbulence to rearrange it- 
self and for new GMCs to collapse with the flow crossing 
time and introduced the dimensionless parameter x ~ 1 to 
para mete rize the uncertainty in the scaling. Using equations 
( W§ and (l3]), we find that t flow = s(Q)tfi BC , where 

*(Q)= ^ ' , (29) 



and thus 



7T ^1/4^/3^' 



,Vtot _ -Mg/coll (on\ 

The collapse fraction / co n depends on the Mach number 
of the turbulence (through the PDF of gas density fluctua- 
tions), the epicyclic frequency of the disc (which determines 
the rotational support on scales ~ h), and the disc Q pa- 
rameter (which quantifies the stability to gravitational col- 
lapse). It also in general depends on the excitation of spiral 
structure and other large-scale disc disturbances by internal 
dynamics or interactions with external perturbers. 



If Q < 1, the disc is Toomre-unstable and / co ii ~ 1- For 
Q > 1, a quiescent disc supported by thermal pressure would 
be stable and / co ii — > 0. However, a realistic galactic disc is 
subject to turbulent fluctuations, so that even if Q > 1 when 
evaluated using averaged properties (e.g., over a ring of finite 
width in the galaxy), random fluctuations imply a certain 
probability that Q < 1 at some locations. Spatially-resolved 
spectroscopy of z ~ 2 star-forming galaxies in fact indicate 
that the locations of star-forming clumps correspond to local 
minima in Q, wher e Q < 1, while Q often exceeds unity 
outside the clumps ( |Genzel et al.^OllI l P 7 ] 

In Appendix[C] we provide a simple heuristic derivation 
for /coii(Q) and summarize the more detailed calculation 
from |Hopkins| ( |2012a[ ). However, these results are sensitive 
to the amplitude of fluctuations on a smoothing scale ~ h 
(identified with GMCs) and are subject to significant uncer- 



10 Since the Ho observations trace star formation, these maps 
likely miss a substantial area where Q > 1. 
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tainties because the properties of the large-scale turbulence 
depend on the mechanism driving the turbulence and the ef- 
fects of disc rotation. Here, we carry out analytic estimates 
using a parameterization of f co ii(Q) that captures its essen- 
tial behavior: 



/coii(Q) = PQ~ 



(Q>1) 



(31) 



The results of Appendix [C] suggest that reasonable param- 
eters are ft ~ 0.5 and a ~ 5 — 3, for Mach numbers 
M = 10 — 100, respectively. 

If GMCs are in steady state in the sense that they pro- 
cess gas at the same rate as they are supplied with gas, 
Mq^ c = Mqm C /toMC- This is likely usually a reasonable 
assumption in practice, because we expect the lifetimes of 
GMCs to be comparable to the free fall time of their host 
disc (4gmc ~ 1), an d the latter is a lower bound for the time 
scale over which the host disc changes. Then 



IGMC 



and thus equation ( 25 1 implies 



GMC 

xs(Q) £int ' 



(32) 



(33) 



Using c T = Q/gff/2 (eq. [9}, equations |22|, pl| and p3 
can be combined to eliminate e| al and solve for Q: 



/3P (P*/m*) GMC 



l/O + l) 



(34) 



_ 2x^/4 /g(T 

(This solution assumes that is a constant; for the case of 
a pure gas disc with no external potential, <f> ~ l/Q (Sj2| 
and the solution involves an additional power of Q.) We can 
use this solution to directly relate the galaxy-averaged and 
GMC star formation efficiencies: 

, (a-l/2)/(«+l) 



4 x 2V 4 x 



1/4 



/go- 



P/3 (P*/m*) 



(35) 

x (/3 £ p n r c ) 3/2(a+1) 



The limit a — > oo, in which the probability of GMC 
formation is suppressed very rapidly as Q exceeds unity, is 
particularly enlightening. Then, 



1 



and 



1/2 



2 5 / 4 J" (P*/m*) (P*/r 



'(p*A 



(36) 



(37) 



In this limit, in which Q approaches exactly unity, the 
galaxy-averaged star formation efficiency per free fall time 
is set by the dimensionless ratio /g(7/(P*/m*), a proxy for 
the ratio of the strength of gravity to the strength of feed- 
back. Perhaps most importantly, in the a — > oo limit, e| al is 
completely independent of how efficiently gas is converted 
into stars once collapsed in GMCs. The rate limiting step 
for star formation in this case is the rate at which GMCs 
form in the galactic disc. For Q = 1, ct = /gf/2 (eq. [9| and 
thus £g al is also proportional to the ratio ct / (P* /m*). 

For finite a > 1, Q is regulated to a value that can 
exceed unity and there is a weak dependence of e| al on e^^ 10 . 
This is illustrated in Figure |2j in which we show equilibrium 



so luti ons for Q , /gmc, and e| al (simultaneously satisfying eq. 



a 



22 



and 



33 1 as a function of e^ t obtained using equation 



( C7 1 for the collapsed fraction as a function of Q instead of 



the power-law approximation in eq. (31 1. The equilibrium 
solutions assume a constant £gmc = t g lsc (£gmc = 1; see 



[2.2|. The left panel in Figure [2] is representative of a low- 
mass or gas-poor galaxy with f g a = 20 km s _1 and the 
panel on the right is illustrative of a massive gas-rich system 
(such as an ULIRG or a high-redshift star-forming galaxy) 
with f g a = 100 km s -1 . In the limit egf c -> 0, GMCs 
do not form any stars and the assumptions of our model 
break down. We thus show the equilibrium solutions only for 
ef^ 1 ^ 0.02, the minimum efficiency implied by equation 
(38 1 below for E g ^ 100 Mq pc -2 . The disc-averaged star 



formation efficiency per free fall time increases slowly with 
the integrated GMC star formation efficiency, while /gmc 
decreases. For example, for a Mach number M — 30, e| al 
increases by a factor ~ 3 when e^^ 10 increases by a factor of 
20 from 0.05 to 1. This is realized by the global Q parameter 
increasing above the classical stability threshold of unity, so 
that GMCs form only where turbulent fluctuations cause the 
self-gravity of the gas to exceed its turbulent and rotational 
support. Figure [2] also shows that for small o^f~ c , the star 
formation rate is too low to support the disc to Q > 1 and 
the disc settles to a Q < 1. 

In Figure [3] we show the same equilibrium solutions but 
as a function of E g for a simple model of how e^ 10 depends 
on Eg for GMC disruption by radiation pressure on dust. 
Specifically, we assume 



GMC 

£int 



. f ttGE k c 



0.35 



(38) 



els of Murray et al 

GMC 



based on the scaling arguments and 1-D numerical mod- 

These authors showed that 
for GMCs that are opti- 



(2010) 



eSZ™ ~ (ttGEgmcc)/ 2(L/M*)] 
cally thin to far-infrared radiation (where Egmc is the GMC 
surface density), reaching a constant ~ 0.35 as GMCs be- 
come optically thick to far-infrared photons. Here, we iden- 
tify Egmc with E g , the mean gas surface density in the 
galaxy, and adopt a light-to-mass ratio L/M t = 3,000 cm 2 
s -3 . In this case, Figure [3] shows that e| al increases by a fac- 
tor of a few with E g ranging from 100 to 1, 600 Mq pc -2 , but 
that /gmc decreases with E g over this range. The flattening 
at Eg = 1, 600 Mq pc -2 is caused by the saturation of e^it 10 
in the optically thick limit. We return to the prediction of 
decreasing /gmc with increasing E g in 94.51 



4 COMPARISON WITH OBSERVATIONS 
4.1 CO conversion factor 

One complication in comparing star formation models to ob- 
servations arises because the molecular gas mass is usually 
estimated using CO emission lines that are optically thick. 
Furthermore, the CO mass is a small fraction of the total 
molecular mass and is therefore only an indirect tracer. The 
Xco = Nu 2 /Ico conversion factor between CO intensity 
and molecular gas column density depends on several fac- 
tors, including gas density, gas metallicity, gas temperature, 
and the ambient radiation field (e.g., Maloney fc Black|1988 



|Genzel et"aT]|2012| |Narayanan et al.||2012| |Feldmann et al.| 
|2012[ ). An equivalent factor aco = Eh 2 //co is related to 
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Xco by Xco = 
cm -2 / (K km s~ 



6.3 x 10 19 Qco, where Xco has units of 
L ) and qco has units of Mq pc _2 /(K km 



In the local Universe, observations find a nearly con- 
stant qco = 3.2 for ordinary galaxies, including the Milky 
Way (e.g.,|Strong fc Mattox||1996| |Dame et alT||2001| |Blitz 
|et al. 12007] ). For local ULIRGs associated with galaxy merg- 
ers, a significantly lower conversion factor aco < 1 is in- 
ferred on the basis that a Milky Way-like factor would im- 
ply a gas mass exceeding the dynamical mass of the galaxy 
I Solomon et al.|1997[ ). On the other hand, qco rises at low 



metallicities such as in the Small Magellanic Cloud (Blitz 
|et al.|[2007[ ) and for high-redshift galaxies below a criti- 



cal stellar mass ( Genzel et al.|[2012 1 , probably because CO 
molecules are photo-dissociated by the UV radiation field 
as the attenuation by dust decreases (e.g., [van Dishoeck fc| 
Black|1986| |Wolfire et al.|2010 |. 

The CO conversion factor is important for our compar- 
ison to observations because different assumptions lead to 
different E* — E g relations and galaxy-averaged star forma- 
tion efficiencies e| al . One common assumption, adopted for 
example in the observational study of |Genzel et al.| ( |2010[ ) , 
is that of a bimodal conversion factor qco = 3.2 for non- 
mergers and qco = 1 for mergers. This assumption leads to 
a E* — Eg relation that is also bimodal, with mergers occu- 
pying an elevated locus. A similar result was found by |DaddT| 
et al. (2010a I. Theoretical considerations have led other au- 



thors to propose an qco factor that depends smoothly on 



galaxy properties. Ostriker & Shetty (20111 suggested that 



Qco should depend primarily on E g . Interpolating between 
ceo = 3.2 (non-merger value) at E g = 100 Mq pc -2 and 
qco = 1 (merger value) at E g = 1,000 Mq pc -2 with a 
power-law, they propose that 



"CO 



(39) 



1,000 M pc- 2 , 
Using this expression, they showed that the merger and non- 



merger galaxies in the Genzel et al. ( 2010 1 sample align on 
an uni-modal E* — E g relation. In this picture, mergers have 
lower qco because they have higher characteristic E g . 



Narayanan et al. (2012) performed radiative transfer 



calculations on simulations of merging and non-merging 
galaxies and predicted the dependence of qco on galaxy 
properties. Based on these calculations, these authors also 
advocate an qco that varies smoothly with galaxy prop- 
erties, rather than a bimodal distribution between mergers 
and non-mergers. For solar metallicity, the calculations of 
Narayanan et al. (20121 yield a dependence on E g close to 



the empirical fit of Ostriker & Shetty (20111 in equation 



(391. Narayanan et al. (20121 also show, in agreement with 



Ostriker & Shetty (2011 1, that assuming a smoothly- varying 
Qco as predicted by their analysis leads to a E* — E g re- 
lation that has a smaller scatter than when assuming a bi- 
modal model. They further argue that the scatter is reduced 
even for a range of galaxies excluding mergers, so that the 
smaller scatter is not simply due to using an unimodal Qco, 
indicating that the dependence on E g is important. 

A recent observational study using dust masses to esti- 
mate the gas masses of z ~ 2 galaxies ( Magdis et al.|[2~012 l 
infers CO conversion factors consistent with the metallicity- 
dependent theoretical prediction of Narayanan et al.| ( f20T2t ) 
within uncertainties. This is the case even though the Qco 
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Figure 4. Disc-averaged star formation law. All data are from 
the compilations of star-forming galaxies of Genzel et al. (2010) 
and Tacconi et al. (2012). We distinguish between local z ~ 
galaxies (blue squares) and high-redshift z ~ 1 — 3 galaxies (red 
triangles). For each redshift interval, the data include both merg- 
ing and non-merging galaxies. The gas surface densities were con- 
verted from the values reported by Genzel et al. (2010) and Tac- 
coni et al. (2012) using an aco factor interpolating smoothly 
between standard merger and non-merger values using equa- 
tion ( | 39 [ t - We assume that all gas is molecular. The solid b lack 
line shows the theoretical prediction S* oc S| in equation a 1 S b 
obtained by balancing the momentum input from stellar feed- 
back with the vertical weight of the disc gas. The line assumes 
Q = <f> = 1, Ft/m* = 3,000 km s _1 , and the dimensionless 
normalization T = 2. The parallel grey lines indicate the range 
Q ~ 0.5 — 1.5 expected for gas-rich galaxies (see Fig. [2J and the 
dashed segments show extrapolations to S g < 100 Mq pc — 2 . 



factor of ordinary high-redshift star forming galaxies is also 
formally consistent with a standard Milky Way value, owing 
to two factors. First, while ordinary high-redshift star form- 
ing galaxies have elevated star formation rates relative to lo- 
cal galaxies, their star formation is typically more spatially- 
extended than that in local mergers, so that the gas surface 
densities are often not as high as those of local mergers. Sec- 
ond, high-redshift galaxies generally have lower metallicity 
Z, which [Narayanan et al. ( 2012[) predict implies a larger 



qco oc 1/Z at fixed E g . 

In light of these considerations, we use an qco depend- 
ing smoothly on E g as in equation ( 39 1 in our analysis. 



4.2 Disc-averaged star formation law 



In Figure)^ we plot E* versus E g for the data compiled by 



Genzel et al. ( 2010 1, supplemented by a larger sample of 2 - 
1 — 3 star- forming galaxies from |Tacconi et al.| ( |2012| . This 
compilation includes both ordinary and merging galaxies. 
In accordance with the range of applicability of our model, 
we focus on data for galaxies with E g > 100 Mq pc -2 and 
assume that all gas is molecular. We convert the gas surface 
densities from those reported by |Genzel et al. (20101 and 
|Tacconi et al.| ( [2012[ ) (who assume a bimodal Qco) to values 
obtained using equation ( 39 1. Figure[4]shows the data points 
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Figure 5. Disc-averaged star formation efficiency per free fall time as a function of isothermal potential velocity dispersion (<r = u c /v2, 
where u c is the maximum circular velocity), gas mass fraction / g , gas surface density E g , and dynamical time in the disc (t^ yn = R\/2/ v c, 
where R1/2 is the half-light radius). The observed disc-averaged star formation efficiency per free fall time is non-universal and increases 
with increasing gas fraction (upper right panel) , as predicted by our feedback model (eq. |37} . The overall scatter in the star formation 
efficiency from galaxy to galaxy is a factor ~ 100. In each panel, the scatter relative to the best-fit linear relation in log-log space is 
quantified by the stddev statistic described in j ]4.3| The data are taken from the observational compilations of Genzel et al. (2010) and 
Tacconi et al. (2012) but converted to a smoothly-varying CO conversion as described in ^4] We distinguish between local z ~ galaxies 
(blue squares) and high-redshift z ~ 1 — 3 galaxies (red triangles) as in Figure [4] In accordance with the assumptions of our modeling, 
we only show data for galaxies with E g > 100 Mq pc -2 . Some estimated gas masses exceed the galaxy dynamical mass (/ g > 1) by a 
factor up to ~ 3, but this does not significantly affect our conclusions (see j]4,3[|. 



converted in this way lie on a well-defined unimodal E* — E g 
relation. 

On this plo t, w e also show the model prediction oc 
E g in equation (18 1. We set Q = = 1 in this comparison 
and also hold P*/m* fixed to the fiducial value of 3, 000 km 
s _1 . Departures from these assumptions are encapsulated 
in the dimensionless factor T; a good match to the data is 
obtained for T = 2. The fact that a good fit is found for 
a constant P+/m+ is consistent with the weak dependence 
of the momentum input by SNe on ambient density and 
turbulent gas velocity dispersion discussed in Appendix [A] 
Interestingly, the star formation law that fits the E g ^ 100 
Mq pc -2 data also fits the data well for E g as low as 30 Mq 
pc -2 , suggesting that our model may apply somewhat more 



broadly. Parallel grey lines indicate the range Q w 0.5 — 1.5 
expected for gas-rich galaxies (see Fig. [2}. 

At very high gas surface densities, supernovae become 
increasingly ineffective due to radiative losses while radi- 
ation pressure on dust becomes more important in the 
volume-filling ISM (Appendix|A|. In the limit of the disc be- 
ing optically thick to reprocessed infrared photons, P*/m* oc 
E g , so we expect a flattening to E* oc E g . This is not vis- 
ible in Figure [4] because of the paucity of the data points 



for Eg > 10 Mq pc (eq. A9 1 but is relevant for galactic 



nuclei ( Thompson et al.|2005[). Observati ons of a sample of 
local active galactic nuclei by |Hicks et al. ( 2009[ ) in fact sup- 
port such a flattening in the E* — E g relation at E g ~ 10 4 

Mq PC" 2 . 
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Figure 6. Disc-averaged star formation efficiency per free fall 
time as a function of gas mass fraction, assuming a bimodal ceco 
conversion factor. The data are the same as in the upper right 
panel of Figure [5] but assume a bi-modal «co f° r mergers and 
non-mergers, as in Genzel et al. (2010) and Tacconi et al. (2012), 
rather than the continuous model in equation ( |39| . This results 
in a larger scatter in the data and a larger number of inferred 
gas masses exceeding the dynamical mass of the galaxy. As in 
the case of the continuously-varying ctco assumption, the data 
do not strongly support the existence of a universal value of e^* 1 . 



4.3 Disc-averaged star formation efficiency 

It is also useful to compare the star formation efficiency per 



free fall time in the disc predicted by equation ( 35 1 to the 



observations. We again compare with the galaxies compiled 
by |Genzel et aL] (|2010[> and |Tacconi et~aI1 (|2012|), but re- 



strict ourselves to systems with E g > 100 Mq pc in the 



rest of this section. Genzel et al. (20101 and Tacconi et al 



(20121 calculate the star formation efficiency relative to the 
dynamical time of the disc, defined as tdyn = R/v c , where 
observationally R is taken to be the half-light radius and v c 
is the maximum circular velocity. Since our star formation 
efficiency is defined relative to the free fall time in the disc 
(201, we multiply the values reported by these authors by a 
factor of 1.14 to use a consiste nt t ime scale (this is the ratio 



t ff lsc /tdyn evaluated using eq. (131, for Q = cj> = 1) 



We compare the observationally inferred e| a against 
several independent variables in order to investigate the rel- 
ative scatter of different relations. In Figure ^ we show 
as a function of the observationally-inferred a, / g , E g , and 
idyn- The isothermal potential velocity dispersion is esti- 
mated as a = v c /\/2, where v c is the maximum observed 
circular velocity. The gas mass fraction is estimated from 
the observations as 



k = 



/2 



2a 2 



(40) 



where R1/2 is the half-light radius. This follows from the 
definition / g = M g (< Ri/2)/M to t{< R1/2) and the simple 
dynamical model in §2.1| This approach is not guaranteed 
to yield an observationally-estimated / g < 1, as must be the 
case physically. Figure [5] shows that it results in a subset 
of galaxies with inferred / g w 1 — 3. We attribute this pri- 



parameters entering in equation ( 40 1 and to the simplified 



dynamical model used. This small subset of galaxies with 
/ g > 1 does not significantly affect our discussion. We em- 
phasize that it is not due to our assumption of an aco fac- 
tor varying continuously with E g . Adopting a bimodal aco 
with ordinary star-forming galaxies at high redshift having a 



Milky Way-like value as in some other studies (e.g., Daddi et 
aT1|2010b| |Genzel et aT|2010| |Tacconi et aL|2012| ) increases 



the inferred gas mass for such galaxies and thus increases 
the number of systems with / g > 1. This is shown in Figure 
[6] which shows the same data as in the upper right panel 
of Figure [5] but assuming a bi-modal aco for mergers and 
non-mergers, as in Genzel et al. (2010) and Tacconi et al. 



(2012), rather than the continuous aco from equation (39 1. 

For each relation in Figure [5] we show in the top left 
corner of the panel the statistic stddev calculated by first 
fitting a linear relation in log-log space, then evaluating the 
standard deviation of the data points from the best fit. The 
uncertainty on the data points (not shown) is dominated by 
systematics; since it is difficult to estimate accurately, we 
assign equal weight to each data point. 

The smallest scatter for the comparison of e| al as a func- 
tion of a single independent variable in Figure [5] is relative 
to / g , for which there is a monotonically increasing relation 
with stddev=0.39. For the other variables (a, E g , and tdyn), 
the scatter is larger stddev=0.48-0.54. Figure [6] also shows 
that the scatter in e| al versus / g is substantially increased 
when a bimodal aco is assumed (stddev=0.55). This is, we 
believe, independent evidence that a continuous aco (as in 
eq. 39 1 is a better approximation. Regardless of whether a 



continuous or bimodal aco is adopted, the inferred values of 
e| al vary by more than a factor of 100 for the molecule-rich 
galaxies in the sample considered. We return to this point 
in So~2l 



In Figure [7] we show similar comparisons but as a func- 
tion of combinations of / g , a, and E g . The left panel shows 



e| al as a function of f g a (stddev=0.36). The solid line in 
this panel is the theoretical prediction e| al oc f g a in equation 
(37 1, valid in the limit a — > 00, for 0=1, P*/m* = 3,000 
km s _1 , and T = 2 (as in Figure [4]). The agreement be- 
tween the prediction and the data (both in normalization 
and in slope) is reasonable given the uncertainties in the 
data and the simplications made in deriving this theoret- 
ical prediction. Moreover, this agreement emphasizes that 
the large scatter in e| al versus E g in Figure [5] is not simply 
scatter about a universal value of e| al . Rather, the scatter 
is significantly reduced when the dependence of £g al on f g a 
is accounted for. For finite a, as must be the case in re- 
ality, the upturn from this simple prediction at large / g cr 
suggested by the data might be explained by an increase of 
^int 10 with E g (see Fig. |3p or a decrease of P+/m* with gas 
density, both of which are expected theoretically. The right 
hand panel of Figure [7] shows e| al as a function of the best- 
fit linear combination of log 10 f g a and log 10 E g . The best fit 
depends weakly on log 10 E g and the stddev=0.34 statistic is 
marginally improved relative to the fit with respect to / g cr 
alone, indicating that f g a is the most important parameter. 
While the simple GMC disruption model in equation 



(38 1 predicts a modest increase of e| al with E g at fixed f g o~ 
(see Fig. [3J, a trend of e| al with E g is not apparent in Figure 
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Figure 7. Same as in Figure [5] but for the disc-averaged star formation efficiency per free fall time as a function of combinations of / g , 
a, and Eg. The solid black line in the left panel shows the theoretical prediction e^ 1 <x f s a of our feedback model in equation II37L valid 
in the limit a — > oo, for </>=l, P*/m* = 3, 000 km s~\ and T = 2 (as in Fi gure 4|. The observed disc-avcraged star formation efficiency 
per free fall time is non-universal and increases with / g <r, as predicted by our feedback model (eq. |37| ), This agreement emphasizes that 
the large scatter in e^ al versus £ g in Figure [5] is not simply scatter about a universal value of , since the scatter is significantly 
reduced when the dependence on f g a is accounted for. The horizontal axis in the right hand panel is the best-fit linear combination of 
log 10 /gcr and log 10 S g . The scatter is not significantly reduced by allowing a dependence of the disc-averaged star formation efficiency 
on S g , indicating that f s a is the most important parameter. 



[5] However, the predicted trend as a function of E g is weak 
relative to the scatter in the data points in Figure [5] More 
importantly, the trend is only predicted at fixed / g <r, while 
the data points cover S3 1.5 dex in f g a. The prediction of 



increasing e|. al with E g also relies on the accuracy of the 



GMC disruption model, which equation (381 undoubtedly 



oversimplifies. The more robust and general prediction of 
our feedback-regulated model of star formation is that e| al 
scales with f g a (eq. 37 1, which the data support (Fig. [7|. 



4.4 Turbulent velocity dispersion 

Our theory also predicts the turbulent velocity dispersion 
ct as a function of galaxy properties. In the limit a — > 00, 
Q — > 1 and equation Q implies that 

2 



Ct 



(0=1). 



(41) 



In the same limit, Ct oc e| a (eq. 221, a result previously 



noted by Ostriker & Shetty (2011 1. The latter authors inter- 
preted this result in the context of a fiducial value e| al ~ 0.01 



motivated by a combination of observations ( |Krumholz fc 
Tan||2007| [Bigiel et al.||2008| |Evans et al.||2009| |Blanc et ~ 



2009 1 and numerical simulations of turbul ent gas (|Krumholz 
fc McKee|2005] |Shetty fc Ostriker|2TjT2| ). [OltrTker fc Shetty 
(20111 pointed out that in their numerical simulations ct 
depended very weakly on E g , in agreement with the lack of 
trend of e| al with E g in the observational data compiled in 
Figure [5] 

Nevertheless, since different galaxies have different / g cr 
and £g a , we expect ct to vary from galaxy to galaxy. This is 
supported by observations. In particular, local U LIR Gs have 
a « 200 km s _1 and / g S3 0.5, so that equation (41 \ implies 



that ct S3 50 km s" 1 in such systems, comparable to what 
is inferred observationally ( |Downes fc Solomon|1998"l ). High- 
redshift (z ~ 2) star-forming galaxies have comparably large 
gas fractions (e.g., Daddi et al. 2010b; Tacconi et al.|[2010[ ) 
and are also inferred to have velocity dispersions that are 

Genzel et al. 2011[) relative t o 
3 10 km s" 1 i Dib et al.|2006| . 



elevated ( Cresci et al. 2009 



local spirals, which have Ct 



4.5 Fraction of disc gas mass in GMCs 

A key aspect of the physical picture presented here is that, 
viewed from the GMC perspective, the star formation rate 
is determined primarily by the rate at which GMCs form in 
galactic discs and is nearly independent of sub-GMC scale 



physics. In particular, provided that a> 1 (see eq. 31 1, the 
galactic star formation efficiency per free fall time e| al de- 
pends only weakly on the efficiency with which GMCs turn 
their gas into stars (thit 10 )- This is possible because for any 
e int I ° t ne GMC formation rate adjusts so as to maintain 
the value of e| al set by the balance between momentum in- 
jection from stellar feedback and the gravitational weight 
of the disc gas at the mid-plane (eq 



18) 



For any model of 
the dependence of the GMC lifetime ^gmc = iGMcift' sc and 
^int 10 on galaxy properties, this translates into a prediction 
for /gmc, the mass fraction of the disc gas collapsed into 
gravitationally-bound GMCs at any given time (see Fig. [2] 
and [3] for specific examples). 

As discussed in ^2j GMCs are defined in our model as 
gravitationally-bound clouds. This definition does not neces- 
sarily coincide with clouds identified based on their molec- 
ular gas content in the Milky Way or other Local Group 
galaxies. The properties of GMCs defined in this way also 
differ from clouds identified in quasi-virial equilibrium. In 
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general, the mass fraction of gas in gravitationally-bound 
GMCs is greater than the mass fraction in quasi-virialized 
clouds, and their lifetime as gravitationally-bound clouds is 
also longer than the lifetime as a quasi-virialized clouds. In 
a purely molecular ISM, gravitational boundedness provides 
a well-defined criterion for identifying GMCs. When obser- 
vations only permit measurement of the mass fraction in 
a more advanced state of gravitational collapse, our predic- 
tions for /gmc should however be interpreted as upper limits 
to the mass fraction in such dense clouds. It is worth noting 
that /gmc (or a variant) is generally an input into analytic 
models (e.g., Krumholz et al. 2009|), but in this work we 



derive /gmc and show how it depends on galaxy properties. 

We consider now the simple scenario described more 
quantitatively in §3.3| and Figure [3] in which the GMC 
lifetime scales with the disc free fall time and GMCs are 
dispersed by radiation pressure on dust. In this scenario, 
e iSt LC saturates at a value ~ 0.35 as GMCs become opti- 
cally thick to far-infrared photons. Quantitatively, this ar- 
gument depends on the detailed properties of GMCs and 
how effectively re-processed far-infrared radiation is trapped 
by scattering on dust grains. Thanks to its scaling with 
Egmc, radiation pressure on dust is the only known stellar 
feedback process potentially capable of disrupting GMCs in 
the densest starbursts such as Arp 220. Thus, if radiation 
pressure on dust were less effective than estimated by |Mur-] 
ray et al.| (2010| ) (e.g., |Krumholz fc Thompson||2012"l ), then 



GMCs in luminous starbursts would likely have even larger 
fint 10 ~ A robust prediction is therefore that in lumi- 
nous starbursts e S^ 10 should exceed the Milky Way value 



Williams & McKee 



1997 



Murray 20111 



eF n f° ~ 0.1 (e, 
by a factor ~ 3 — 10. 

Consider the example of the local ULIRG Arp 220 
(E g ~ 10 4 M ; Scoville et al. 1997). We fiducially assume 
£gmc ~ 1 and estimate e| al « 0.05 using the data compiled 
. Genzel et al. fl2010|). For th is galaxy, the GMC disruption 



theory of |Murray et al.| ( |2010[ > predicts e^f » 0.35. There- 
fore, equation (25 I implies /gmc ~ 0.14, in agreement with 



the more accurate numerical solutions in Figure [3] 

Interestingly, high-resolution aperture synthesis CO ob- 
servations of the nuclear gas disc in Arp 220 indicate a high 
area filling factor ~ 0.25, suggesting that the molecular gas 
is more uniformly distributed than in less extreme systems 
( Scoville et al.|[l997 l. In contrast, the volume-filling factor 
of molecular clouds in the inner Galaxy is ~ 0.005 (Roman- 



Duval et al. 2010 Kennicutt & Evans 2012 1 and is much 



lower in the rest of the Milky Way. The apparently-smooth 
CO-emitting gas distribution in ULIRGs in general has been 
interpreted to be connected to the smaller aco conversion 
factor in those systems (e.g., Solomon et al.|[l997 Downes 



fc Solomon|1998 1, in qualitative agreement with the predic- 



tion that /gmc is relatively small. It must be noted, though, 
that there is a possible alternative interpretation. In the cen- 
tral ~ 100 pc of ULIRGs, the mean gas density is typically 



11 Powerful outflows from star- forming clumps have been de- 
tected in high-redshift star-forming galaxies | |Genzel et al.|2011| 
|Newman et al.|2012| , suggesting that stellar feedback is effective 
at dispersing them. Winds driven by young super star clusters 
are also observed to unbind cluster gas in the prototypical merg- 
ing g alaxies NGC 4038/4039 (the Antennae; |Gilbert fc Graham| 
20071. 



nu ~ 10 4 cm -3 (see eq. ^ and nearly all the gas is in molec- 
ular form. Thus, even if a dominant mass fraction of the gas 
is collapsed in gravitationally-bound GMCs (/gmc ~ 1), it is 
conceivable that the CO-emitting area filling factor appears 
unusually large simply because of the inter-cloud medium is 
sufficiently dense to be optically thick. More detailed mod- 
eling is needed to distinguish these possibilities observation- 
ally. 

High-redshift (z ~ 2) star-forming galaxies, with circu- 
lar velocities and gas fractions comparable to local ULIRGs, 
also appear to have relatively low /gmc- Generally, /gmc is 
not directly measured since most existing observations in- 
stead trace the star formation rate or stellar mass of high- 
redshift galaxies. Estimates based on cosmological simula- 
tions suggest /gmc ~ 0.1 — 0.2 (ICeverino et al. 20101, al- 



though in our picture this quantity is sensitive to feedback 
parameters, which are highly uncertain in existing simula- 
tions. Recent observations with the Hubble Space Telescope 
however indicate that at most a fraction ~ 0.02 — 0.07 of the 
total stellar mass of such galaxies is contained in clumps. 
Integral field measurements by Genzel et al. (2011 1 further- 



more find no prominent kinematical imprint at the locations 
of star-forming clumps, indicating that they cannot be dom- 
inant by mass. 

The nuclei of local gas-rich spirals offer another promis- 
ing opportunity to test our /gmc predictions. Galactic nuclei 
can have gas surface densities E g > 100 Mq pc -2 and have 
the advantage of being numerous in the nearby Universe. 
They can thus be observed at spatial resolution sufficient to 
resolve their molecular clouds. Rosolowsky fc Blitz| ( |2005[) 
observed the nucleus of M64, where S g w 160 M© pc~ 
over a ~ 2 kpc diameter region. In spatially-resolved studies 
like this one, optically thick emission from 12 CO can be used 
to infer to total gas mass, while transitions of molecules such 
as 13 CO (only marginally optically thick) and HCN (with a 



critical density for the J 



transition of riH 



10 s 



cm" ) can be used to identify the locations of GMCs. The 
masses of GMCs can then be derived from measurements of 
their velocity dispersions and radii. In the future, such stud- 
ies could be extended to statistical samples using existing 
interferometers such as the Combined Array for Research in 
Millimeter- wave Astronomy (CARMA^and ALMA. 

Independent support for our physical picture is pro- 
vided by numerical simulations. |Dobbs et al. (20111 per- 
formed simulations of galactic discs with supernova feed- 
back and varied the feedback energy per molecular gas mass 
formed via their e parameter. In an experiment in which 
e was increased from 0.05 to 0.2, the fraction of disc gas 
mass contained in GMCs is reduced by a factor of approxi- 
mately two. For a fixed gas mass in which the feedback ki- 
netic energy is injected, the momentum returned per GMC 
mass formed scales as oc y/e in the simulations of Dobbs 



et al. (2011 1. Thus, the experiment of Dobbs et al. (2011 1 is 



consistent with the prediction in equation ( 25 1 that /gmc 
should be reduced oc (e£f c ) _1 . 



http:/ /www. mmarray.org/ 
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5 DISCUSSION 

5.1 Summary of the main results 

We presented a physical picture connecting star formation 
in giant molecular clouds with the global star formation rate 
in galaxies. We focused on galaxies with gas surface density 
100 Mq pc~ 2 to simplify the theoretical treatment 
by avoiding explicit consideration of a multiphase (atomic 
and molecular) ISM; we expect, however, that our quali- 
tative points apply to galaxies with lower surface density 
as well. The high gas surface density regime includes local 
starbursts, most star-forming galaxies at redshift z > 1, and 
the nuclei of nearby ordinary galaxies. Our theory builds on 
previous models of star formation which focused on feed- 
back on scales of galactic discs (e.g., [Silk 1997; Th ompson 
et al.|2005||Ostriker fc Shetty|20iT} or GMCs (e.g.,|Matzner 



2002 Krumholz fc Matzner 2009 Murray et al. 20101. In 



particular, our model explicitly determines the relationship 
between the formation and destruction of GMCs and the 
galaxy-averaged star formation rate. 

In our theory, the star formation rate of a galaxy and 
thus the star formation efficiency per free fall time e^ 3,1 is 
set largely by the balance between the gravity acting on the 
disc gas and the strength of stellar feedback in the bulk of 
the ISM (see also Thompson et al. 2005; Ostriker & Shetty 
2011). If supernovae dominate the turbulence in the volume- 



filling ISM, this yields a star formation law oc S g (eq. 18 1 



The integrated efficiency with which GMCs convert their gas 
into stars, e^it 10 ! a l so enters into the galaxy-averaged star 
formation rate but only weakly (in spite of the fact that 
we assume that all star formation occurs in GMCs). This is 
possible because the formation rate of GMCs adjusts itself 
to ensure that a E* — E g relation set by vertical hydrostatic 
equilibrium in the disc is satisfied, for any reasonable sce- 
nario of sub-GMC physics. This is realized via the Toomre Q 
parameter, which can obtain values above the classical sta- 
bility threshold of unity so that the GMC formation rate is 
regulated. Our physical model is consistent with numerical 
simulations that find that the star formation rate of galaxies 
is not sensitive to the small-scale star formation law (e.g., 
Hopkins et al.|2011 1. 

A key implication of our results is that in very dense 
galaxies where stellar feedback in GMCs is relatively in- 
effective and high efficiencies of star formation in GMCs , 
eSf C ~ 0.35 - 1, are expected (e.g., |Murray et al.|[201o| , 
the star formation efficiency of GMCs can exceed the galaxy- 
averaged star formation efficiency e| al by a substantial fac- 
tor. As a result, the low observed star formation efficiency 
of galaxies does not depend fundamentally on star formation 
being slow within GMCs (although this can be the case as 
well in some galaxies). 

Our theory makes several predictions. In the limit in 
which the GMC formation rate is rapidly suppressed as the 
Toomre Q parameter of the disc exceeds unity (a> 1; see 



eq. 31 1, they are: 



(i) The galaxy-averaged star formation efficiency per free 
fall time e| al oc / g cr/(P*/m*), as implied by hydrostatic 
equilibrium. 

(ii) The turbulent gas velocity dispersion ct ~ /gO"/2 as 
a result of the regulation to Q ~ 1. 

(iii) The fraction of the disc gas mass collapsed in 



gravitationally-bound GMCs /gmc oc (e^t 10 ) -1 m order to 
satisfy the star formation law set by hydrostatic equilibrium 
in the disc. 



In §4.3| we provide strong observational evidence that 
e| al oc / g (T, in agreement with the first prediction (Fig. [7]). 
The second prediction regarding the velocity dispersion de- 
pends only on Q w 1 and is thus very robust. It explains why 
gas-rich systems like local ULIRGs and high-redshift star- 
forming galaxies have velocity dispersions ~ 50 — 10 km 
Downes fc Solomon|l998) |Genzel et al.|201l[ ), sig- 



(e,< 



nificantly above more gas-poor local galaxies like the Milky 
Way (e.g., Dib et al. 2006; G3. The third prediction of a 



relatively small /gmc in dense galaxies where larger GMC 
star formation efficiencies are expected theoretically tenta- 
tively explains the smooth gas distributions and low CO 
conversion factors inferred in local ULIRGs (e.g., Scoville et 
al. 1997, Downes & Solomon 1998) and the small mass frac- 
tion in giant clumps in high-redshift galaxies (Genzel et al. 
2011; Wuyts et al. 2012; p"5L 

We note that our main results can be generalized to 
feedback mechanisms other than supernovae and radiation 
pressure (such as photoionization or cosmic rays; e.g. McKee 
1989, Socrates et al. 2008) by appropriate choices for P*/m* 
and egf c . 



5.2 Role of turbulence 

Turbulence has been invoked to regulate star formation in 
many models. In our theory, galactic discs are supported ver- 
tically by supersonic turbulence and the formation of GMCs 
relies on turbulent density fluctuations. The role of turbu- 
lence in our model is, however, different than in other popu- 
lar theories of star formation. In particular, our theory does 
not rely on supersonic turbulence being maintained within 
GMCs, an assumption that |Krumholz fc McK cc (2005) used 
to predict a nearly universal volumetric star formation law 
of ~ 0.01 of the gas mass turned into stars per free fall 
time in molecular gas (see also |Padoan fc N ordlund| poTT] 
Federrath fc Klessen|2012 1. As mentioned above, our theory 
allows for star formation efficiencies within GMCs as high as 
e inf C ~ £ff MC ~ 1. There is in fact some tentative evidence 
in the Milky Way that some GMCs have star formation effi- 
ciencies significantly larger than predicted by theories based 
on the small-scale properties of supersonic turbulence (e.g., 
Murray|2011 and references therein). We suspect that this 
is even more likely to be true in denser galaxies. 

Even if supersonic turbulence is maintained within 
GMCs, simulations show that in the absence of a mecha- 
nism to support the galactic disc on large scales (such as 
stellar feedback) , self-gravitating regions collapse until their 
free fall time is very short. Therefore, even if the small-scale 
star formation efficiency is low, the star formation efficiency 
per galactic free fall time e| al is ~ 1 in the absence of global 
support (e.g., Bournaud et al.|2010 Hopkins et al.|2012 |. 

Our theory predicts that the volumetric star formation 
law is not universal either on galaxy scales (e| al ) or GMC 
scales (eff MC ). This is seemingly in conflict with the observa- 
tional evidence compiled by Krumholz et al. (20121 in favor 



of a volumetric star formation law with a universal efficiency 
of ~ 0.01 per free fall time. We note that the overall scatter 



in the £g al values inferred by Krumholz et al. (2012 I for the 
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molecule-rich galaxies considered in this work is comparable 
to the factor ~ 100 we find (Figs 5)7 1, regardless of whether 
we assume a continuously-varying or bimodal aco conver- 
sion factor (Fig.|6|. The observational analysis of Krumholz 



et al. (20121 also includes ordinary local group galaxies and 



individual molecular clouds within the Milky Way, and thus 
covers a larger dynamic range than ours. The linear relation- 
ship between E + and E g /ta over this large dynamic range 
likely reflects the fact that the star formation rate on average 
does scale inversely with the local free fall time, as advocated 
by those authors. Averaged over galaxies, however, Figure 
M shows that e| al is not well described by a constant value 
with large scatter. Instead, the star formation efficiency e| al 
correlates well with f g a (particularly for our assumption of a 
continuously- varying aco), in agreement with our feedback- 
regulated theory (eq. |37[ ) but inconsistent with a universal 
volumetric star formation efficiency. 

Finally, we note that our calculations suggest that stel- 
lar feedback can drive the turbulence required to explain 
the observed star formation law and gas velocity dispersions, 
without the need for additional sources such as cosmological 
accretion and radial inflow within galactic discs. Significant 
contributions from such sources cannot, however, be ruled 
out by our arguments. 



5.3 Opportunities for numerical modeling 

The analytic arguments presented here are subject to some 
uncertainties, which could be addressed by numerical model- 
ing. One important uncertainty is exactly how the GMC for- 
mation rate is regulated in a turbulent disc, especially when 
Q > 1 in the disc-averaged sense. We parameterized this 
uncertainty by the dimensionless function / co n oc Q~ a (eq. 
31 1. Analytic estimates suggest that a ~ 2 — 5 ( Appendix [c] 
Hopkins 2012a). However, this is sensitive to density fluctu- 
ations on a scale comparable to the disc scale height. These 
fluctuations in turn depend on how turbulence is driven and 
the effects of disc rotation. 

Similarly, we parameterized stellar feedback processes 
by the effective momentum injected into the ISM per stellar 
mass formed, P*/m*, and carried out our numerical calcu- 
lations using simple estimates for supernovae and radiation 
pressure on dust. These remain uncertain at the factor of 
few level. An important uncertainty for radiation pressure is 
the role radiation-hydrodynamic instabilities in facilitating 
the leakage of scattering photons. This problem has received 
some attention recently (e.g., |Krumholz fc Thompson|2012"| ), 
but direct simulations have only been performed in specific 
idealized settings, such as a laminar background disc. 

While the dynamics of individual supernova remnants 
have been the subject of many studies (e .g.,|Cioffi et al.|1988| 
Draine fc Woods||1991| |Thornton et al.|[l998[ ) the collective 
impact of multiple SNRs in realistic galactic discs does not 
follow simply from these studies. The effective P*/m* de- 
pends not only the ambient density, but also the stage of 
SNR evolution at which it is evaluated. It is not clear how 
to accurately calculate this quantity in a real galaxy with an 
inhomogeneous ISM. Fortunately, the exact value of P*/m+ 
appears to depend only weakly on these details (see Ap- 
pendix El. On the other hand, we also showed in |3.1| that 
the turbulent pressure in the ISM depends on the spatial 



scale on which the turbulence is driven relative to the disc 
scale height, the dimensionless factor fh — L/h. Thus, the 
turbulent pressure depends on how SNRs merge with the 
ambient ISM and/or interact with one another, a problem 
which controlled numerical experiments could address. 
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APPENDIX A: FEEDBACK PARAMETERS 

An important parameter in the calculations of the main text 
is the effective momentum returned to the ISM per stellar 
mass formed, P*/m*. Here we summarize relevant results 
on the importance of supernovae and radiation pressure on 
dust. Other processes - including HII regions, stellar winds, 
and proto-stellar winds - can be important in dwarf or Milky 
Way-like galaxies but are not effective at the densities con- 
sidered in this work ( |Matzner|[2002| |Murray et al.|20"l0 l. 

We distinguish between stellar feedback in the volume- 
filling ISM and in GMCs because the processes that domi- 
nate turbulence driving in the volume-filling ISM are gener- 
ally different from those responsible for disrupting GMCs. In 
what follows, the superscript 'disc' corresponds to a process 
operating directly in the volume-filling ISM, while the su- 
perscript 'GMC corresponds to a process operating within 
a GMC. 



Al Supernovae 

The first SNe explode a time £sN,ist ~ 3.6 Myr following a 
star formation event, which in dense starbursts is after the 



parent GMC has been dispersed (tsN,ist > ^gmc; see eq. 13 1. 
It follows that SNe typically explode in the volume-filling 
phase of the ISM. We are also interested in high-redshift, 
gas-rich star-forming galaxies, for which the lifetime of the 
massive star-forming clumps are inferred to be ~ 100 — 200 
Myr ( |Wuyts et al.|2012| . In these systems, it is not obvious 
that most SNe explode outside their parent GMC. In the 
Milky Way, where GMC lifetimes also exceed tsN.ist, it is 
nonetheless observed that GMCs often start to be disrupted 



prior to the onset of SNe (Rahman & Murray 2010[ ). By 



© 0000 RAS, MNRAS 000, 000-000 



Feedback-regulated star formation 17 



analogy to the Milky Way, we assume that this is also the 
case in high-redshift galaxies. 

In the Sedov and pressure-driven snow plow phases of 
supernova remnants (SNRs), the momentum of the swept 
up ambient medium can reach values exceeding that of the 
original supernova ejecta by a factor 

-1/14 

■>U I — r ! (Al) 



Psr* 



PsN,0 



( nH ) 

VI cm- 3 / 



-1/7 



VSN 



10, 000 km s- 



I Cioffi et al. 1988 1, where _Bsn is the kinetic energy of the su- 



pernova ejecta, wsn is their velocity, 71h is the density of the 
ambient medium, and the exact pre-factor depends weakly 
on metallicity. Noting that the total rate of kine tic energy 
injection by supernovae is Esn ~ 0.01L* (e.g., 



Leitherer 



|et al.|1999| ), where L+ is the stellar bolometric luminosity. 



Psn « 30- 



10 51 erg s - 



-1/14 



( nH ) 

VI cm- 3 / 



-1/7 



(A2) 



We define m* to be the total mass of stars formed per 
supernova event and fiducially adopt m* = 100 Mq (e.g., 
Ostriker fc Shetty||2011[ ). In this case, equation (|A2l) corre- 
sponds to 



■■ 4,800 km s _1 £ 



( nH ) 
V 1 cm" 3 / 



-1/7 



-Esn 
10 51 erg 

VSN 



13/14 



(A3) 



10, 000 km s- 



rdisc 
JSN ! 



where /sn c is the fraction of SNe exploding in the volume- 
filling ISM. We have introduced the dimensionless pre-factor 
£ ~ 1 to parameterize the uncertainty in the normalization 
of equation &A31. In an inhomogeneous medium, the relevant 



density at which to evaluate equation ( A3 1 is the typical 



density where SNe explode. When massive stars outlive their 
parent GMCs, we can assume that SNe explode in random 
locations in the ISM. For a supersonically turbulent medium 
with Mach number M, we define the effective density p e s 
such that 50% of the volume has density < p c g . In Appendix 
[b| we show that p c g w 0.06p(A1/30) _1 ' 2 for a lognormal 
density probability distribution function (PDF). Thus, the 
effectiveness of supernova feedback in a dense ISM is not as 



reduced as would be inferred by evaluating equation ( A3 1 at 
the mean density. 

Equations ( |AHA3[ ) assume a value for the momentum 
boost valid at infinity. Cioffi et al. ( 1988 1 show, however, 



that convergence to this asymptotic value is slow and that 
it is likely not fully realized in practice. In §3.1| we argued 
that for our purposes P*/m* should be evaluated when the 
SNR velocity equals the turbulent gas velocity in the ISM, 
wsnr = ct- We also evaluated P*/m* using the semi-analytic 
approximations of Draine & Woods ( 1991 1 and found that 
P*/m* ranges from 1,500 km s~ to 3,700 km s _1 for 
fin = 100 - 10 4 cm" 3 and final « S nr = 10-50 km s" 1 . The 
density dependence of (P^/m*)^ 1 ^ is thus sufficiently weak 
that to a fair approximation it can be assumed constant for 
most observed galaxies, justifying our fiducial assumption 
of (P*/to*)sn c « 3,000 km s~\ Uncertainty in the precise 
value is encapsulated in the dimensionless parameter T de- 
fined in equation \19\ . 

Using the numerical simulation results of |Thornton| 



et al.| ( |1998[ ), |Thompson et al. ( |2005[ ) estimated P*/m* » 
1,500 km s" 1 (n H /l crn -3 )- 1 ^ 5 for SNe (their eq. 11 and 
associated text, in our notation). This scaling with density 
suggested that supernovae become subdominant relative to 
radiation pressure in the volume-filling dense ISM of star- 
bursts for lower E g than estimated in §A2.2| The |Thornton| 
et al. ( 1998 I results for the late-time momentum of super- 
nova remnants are not, however, consistent with our semi- 



analytic results based on Cioffi et al. ( 1988 1 and Draine & 



Woods (19911. As described above, the latter give a some- 



what larger late-time momentum and a weaker dependence 



on ambient density. In particular, the Thornton et al. ( 1998 \ 



results for the velocity of the supernova shock (their eq. 23, 
which is evaluated at a time of 13to, where to is the time 
of peak remnant luminosity) do not appear consistent with 
those of ICioffi et al.l Jl988| and IDraine & Woodsl (11991 \ 



while other remnant quantities are (e.g., radius, the time to, 
etc.). The origin of this difference is not clear. 

A2 Radiation pressure on dust 

Radiation pressure on dust is the only known mechanism 
capable of disrupting GMCs in the densest starbursts. This 
is because radiation pressure on dust scales with the surface 
density of GMCs, so that it remains important to arbitrarily 



high densities ( |Murray et al.||2010| >. After the CMC is dis- 
rupted, radiation pressure from the surviving massive stars 
acts on the volume-filling ISM. 



A 2.1 GMC disruption 

When evaluating the effective P*/m* from GMC disruption, 
we must account for the fraction of the momentum input 
used to overcome the self-gravity of the GMC. As before, we 
define e^f c as the fraction of a GMC mass that is converted 
into stars and let vt be the terminal velocity of the dispersed 
GMC gas (after being unbound but before being decelerated 
by mass loading of diffuse gas) . Since the dispersed GMC has 
asymptotic momentum Mgmc(1 



GMC 



GMC\ 



-Vt. 



(A4) 



The integrated GMC efficiency e,^ and the terminal ve- 
locity t>{ depend on the properties of the GMC and of the 
emb edded stellar cluster s . 

Murray et al. (20101 report e^^ 10 and V{ values for their 



1-D models of GMC disruption by radiation pressure on 



M 



PC 



dust. For the local ULIRG Arp 220 (E g ~ 10 4 
the local starburst M82 (E g ~ 500 M pc~ 2 ), and the 
z ~ 2 star-forming galaxy BX482 (Eg ~ 200 M pc" 2 ) 
equation {A4} implies (P*/m*)? a d C = 82, 32, and 135 km 
s _1 , respectively. Assuming the scaling in equation I A3 1 
for (P*/m*)iif , (P*/m,)i^ > (P*/m,) r G ^ C unless n H > 
10 11 cm 3 , indicating that supernovae in the volume-filing 
medium dominate over GMC disruption in driving turbu- 
lence in the volume-filling ISM in essentially all realistic con- 
ditions, in agreement with the analytic argument given by 
|Ostriker fc Shettyj poTTj ) . 

Some caveats are in order here, since the effective 
P*/m* from GMC disruption depends somewhat on the in- 
ternal structure of GMCs and their time-dependent evolu- 
tion. We cannot rule out that GMC disruption may be more 
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important for turbulence in the volume-filling medium than 
estimated above. In addition, once the GMC starts to dis- 
rupt, the massive stars continue to input momentum into 
the ambient medium, as discussed in j ]A2.2| This may be 
more important than GMC disruption in driving turbulence 
in the ISmE] 

Giant clumps in high-redshift star-forming galaxies oc- 
cupy an interesting region of parameter space in which 
photoionization is ineffective (Matzner; 2002 Murray et al. 
|2010[ ), radiation pressure is only margi nally capable of dis - 
rupting the clumps (with egf c ~ 0.35; |Murray et al.|2oio| , 
but the free fall time in the clumps 



« 16 Myr 



Afd 



10 9 Mr; 



-1/2 



Be 



1 kpc 



3/2 



(A5) 



is sufficiently long that supernovae may play an important 
role in disrupting them. This should be studied in more de- 
tail in future work. 



A2.2 Direct action on the volume- filling ISM 

The volume-filling medium can also be optically thick to 



reprocessed far infrared (FIR) radiation (e.g.,|Sc ovillc 2003 



Thompson et al.|2005[ |Andrews fc Thompson|2011| )| 14 |The 



momentum flux provided by radiation from massive stars is 
then 



Prad = (1 + Tm) — 



(A6) 



where rnf is the effective IR optical depth. The Tnf term 
enters because confinement of IR photons by multiple scat- 
terings boosts the momentum flux, in a manner analogous 
to the hot gas confinement that gives rise to the momentum 
boost in SNRs. 

The IR opacity of dust is temperature- dependent, but 
peaks at sir ~ 5 cm 2 g" 1 in the model of Semenov et al. 
( [20031 ) for T ~ 100 - i- 000 K («iR dr °P s sharply for T > 
1,000 K owing to dust sublimation). Defining the vertical 
optical depth of the disc as tir, = EgKiR/2 for the case of a 



homogeneous gas distribution, 

2 



TIR 



KlRO_ 

2nGr 
15.4 



(A7) 



KIR 



5 cm 2 g _i y V-^© 

X (2OO km s- 1 



V 100 pc 



where we scale the result with gas-phase metallicity, Z, rel- 
ative to solar. The effective rfpf differs from this estimate 
because anisotropics in the column density distribution can 
provide paths of least resistance through which the scat- 
tered photons can escape. The magnitude of this effect de- 
pends on the column density distribution, and also possibly 
on radiation- hydrodynamic instabilities (e.g., Krumholz & 
Thompson 2012 Jiang et al. 2012 1. We parameterize this 
effect using the dimensionless factor f e g < 1, such that 

T IR = JegTlR. 

Let e ra d ~ 5 x 10 -4 be the efficiency with which star 
formation converts rest mass into radiation. Neglecting con- 
finement of the reprocessed radiation, the radiative momen- 
tum released is simply (-P*/m*) ra d = trade. In general, we 
thus have 



(A8) 



/ 1 1 eff \ _fdisc 
= (1 + Tm )/ rad tradC 

= 150 km s- 1 (1 + ( 

V 5 x 10 



where f^ c is the fraction of the radiation emitted in the 



disc. Equation (A8l equals (P*/m*) S N C when 



10 4 M© pc" 



/eg 

05 



^rad 



5 x 10- 



(A9) 



/■disc 

7sn 

/disc 
rad 



\disc 
JSN 



3, 000 km s" 



Above this gas surface density (corresponding to tir ~ 10 
for the fiducial parameters above), radiation pressure on 
dust dominates the vertical pressure support over super- 
nova feedback. For this optical depth, the analytic estimates 



of Hopkins et al. (20111 imply a modest inefficiency factor 



/eg ~ 0.5. In particular, equation (A9l implies that radia- 



tion pressure is particularly important in galactic nuclei for 
r < 100 pc (as in, e.g., the Thompson et al. 2005 model). 



J In their numerical simulations, jHopkins et al.| ( |2011^ found 
that radiation pressure alone was capable of regulating the star 
formation law to be comparable to that found in Kennicutt's 
(1989) observational sample. We note, however, that the version 
of the KS law |Hopkins et al.| ( |2011} compared to — from |Kenni-| 
|cutt| |1989[ | — is higher in normalization at E g = 100 g cm z by 
a factor of ~ 3 than the data utilized in this paper. As a result, 
the |Hopkins et al.| ( |2011| l radiation pressure feedback only sim- 
ulations of high star formation rate galaxies (e.g., their Sbc and 
HiZ models) would have star formation rates larger than the data 
utilized in this paper by a factor of ~ 3 — 5. However, the subse- 
quent models by Hopkins et al. (2012) include additional feedback 
processes, including SNe, and are in much better agreement with 
the data for the KS law used here. 

14 The ISM is always optically thick in the UV for the galaxies 
we consider. 



APPENDIX B: DENSITY PDF FOR 
SUPERSONIC TURBULENCE 

We assume that the gas density distribution in the ISM is 
determined by the properties of isothermal supersonic tur- 
bulence. Simulations of driven turbulence then imply that 
the mass-weighted density PDF is well described by a log- 
normal distribution: 



fm{y)dy = 



1 



exp 



(y - m) 2 



2ol 



dy, 



(Bl) 



where y = In (p/p) and p is the mean of the distribution 
(e.g.,|Nordlund L Padoan|1999||Ostriker et al.|2001||Lemas-| 
ter & Stone 20081. The constraint that (p) = p implies that 
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Figure Bl. Mass (solid) and volume (dashed) fractions of the 
gas density distribution for driven isothermal supersonic turbu- 
lence for different Mach numbers. These values are calculated as- 
suming a lognormal PDF with parameters determined using the 
prescription of Federrath et al. (2008), for a parameter b = 1/2 
corresponding to a random mixture of solenoidal and compressive 
driving modes. 



p = crg/2. The volume- weighted PDF, fv{y), is given by the 
same expression with the simple replacement p — > — p. 



Federrath et al. ( 2008 1 showed that the lognormal pa- 



rameters depend on both the Mach number M of the tur- 
bulence on the driving scale and on the relative importance 
of solenoidal and compressive driving modes. To a good ap- 
proximation, 



2 
<?0 



ln(l + 6 2 .M 2 ), (B2) 

where 6=1/3 for pure solenoidal driving, 6=1 for pure 
compressive driving, and intermediate values of b correspond 
to solenoidal and compressive mixtures. For a random mix- 
ture of solenoidal and compressive modes, b = 1/2. 

Figure |B1| shows the integrated fractions of the mass 
and volume contributed by over-densities < p/p, as a func- 
tion of p/p. We heuristically postulate that SNe acting out- 
side GMCs are embedded in an ISM of effective density 
[p/p] c a such that Fv(< [p/p]e«) = 0.5. A simple approxi- 
mation is that 



(B3) 



For a starburst with Mach number M = 30, the effective 
volume- filling density is therefore about 6% of the mean ISM 
density. 

APPENDIX C: COLLAPSED FRACTION IN 
TURBULENT DISC 

We derive here an approximation to the fraction of the gas 
mass in the turbulent disc that is unstable to gravitational 



collapse at any time, / co n (see eq. 281. The approximation 



is based on the assumption that most of the collapsed mass 
is contained in the most massive GMCs, corresponding to a 
spatial scale of order the disc scale height h. 

Since /gmc ~ /colliGMC in a steady-state disc with Q ~ 
1 (eq. |32[ | , / co n ~ /gmc if and only if GMCs survive for a 
disc free fall time (£gmc ~ !)• Physically, this arises because 



if GMCs live longer, then the steady-state mass in GMCs 
can exceed the instantaneous fraction of the turbulent disc 
gas mass that is unstable to gravitational collapse. 

For a disc with global Q(disc) > 1, we define the col- 
lapsed fraction / co ii as the fraction of the total gas mass 
with Q(h) < 1, where Q(h) is the Q parameter evalu- 
ated after smoothing the gas density distribution on a scale 
h. In the main text, Q(disc) is abbreviated to simply Q. 
Since there are in general fluctuations on spatial scales > h, 
Q(disc) 7^ Q(h). A more accurate derivation improving the 
standard Toomre analysis for a turbulent disc is possible us- 
ing the excursion set formalism (Hopkins 2012a). For our 
numerical calculations in the main text, we in fact use the 



analytic approximations derived in Hopkins (2012a I, but 



present here a simpler derivation that makes the origin of 
the qualitative behavior with Q(disc) more transparent. 
We must first quantify the gas density fluctuations 



smoothed over the scale h. Equation ( B2 1 was derived for 



simulated turbulence in a static box with periodic bound- 
ary conditions and applies to the point PDF. In general, it 
is non-trivial to generalize this expression to the PDF of the 
density smoothed on a scale ~ h because the latter depends 
on how the turbulence is driven, and how it is suppressed 
on large scales by finite mass and rotation effects. We adopt 
a simple model proposed by |Hopkins| ( |2012b[ ) , in which an 
analog of equation ( B2 1 is assumed to apply on a fc-by-fc 
basis: 



2 



In 1 + b 



M' 2 (k) 



l + 2M 2 (k)/\kh[ 



(CI) 



where M(k) is the Mach number on scale k ~ 1/7?. If the 
turbulence power spectrum E(k) oc k~ p , where c^(k) ~ 
kE(k), then M(k) oc k {1 ~ p)/2 . For a Burgers power spec- 
trum, p — 2, M(k) oc fc~ ' 5 . The real-space variance on 
scale R is then 



2 



dk 2 
k 



(C2) 



-space. 



where we assume a top-hat window function in k 

The denominator 1 + 2M 2 (k)/\kh\ 2 in equation (CI I 



was introduced to capture the suppression of turbulent fluc- 
tuations on scales ~ h owing to a combination of finite mass 
and rotation effects. It is motivated by the Toomre disper- 
sion relation 



2 , 2 
C T K 



27rGE g |fe|, 



(C3) 



which implies that differential rotation in the disc provides 
an effective pressure K 2 k~ 2 analogous to on scales ~ h, 
so that M 2 = Cx/c 2 — > Cx/(c 2 + K 2 k~ 2 ). The suppression 
of the density fluctuations on scales ~ h is critical to the 
strong suppression of / co n for Q(disc) > 1 that we find be- 
low. However, the exact magnitude and functional form of 
this suppression is uncertain and an important problem for 
future work is to quantify this more accurately using numer- 
ical simulations. 

The local Q parameter smoothed on a scale h is 



Q(h) = 



<tct(/i) 
irGp(h)hr 



(C4) 
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Q(disc) Q(disc) 

Figure CI. Left: Fraction of the disc gas mass in a turbulent disc that is unstable to gravitational collapse as a function of the disc- 
averaged Toomre Q parameter. The (overlapping) dashed curves show the heuristic approximation derived in Appendix [c] (eq. |C6| l and 
the solid curves show more accurate results obtained using the excursion set model of Hopkins (2012a) (our cq. |C7| l. These curves assume 
dimensionless parameters a = 1.5 and b = 0.5 in equation |C7| Right: Power-law approximations used for the analytic estimate in §3.3| 
These reproduce the more detailed curves in the left panel reasonably well. 



Taking h as a disc constant]^] Q(h) < 1 if and only if 

J£K > Q(disc). (C5) 
p(disc) 

Thus, the collapsed fraction estimated using this heuristic 
Toomre-scale argument is 

fmiVh > lnQ(disc)) (C6) 



/Toom 
coll 



= -eric 



In Q (disc) — fih 



where the subscript h is used to denote quantities smoothed 
over a scale h. 



Using the excursion set (EPS) formalism, Hopkins 



( 2012a I shows that the "maximum instability scale" in a tur- 



bulent disc is close to but distinct from h. Hopkins (2012a) 



also obtains a more rigorous estimate of the collapse frac- 
tion valid for Q(disc) > 1 and accounting for the full mass 
spectrum of GMCs: 



rEPS 
/coll 



M 

2V1 + M 2 



= erfc 



a In [3Q(disc)/V2] 

■^2 In [1 + 0.5b 2 M 2 /^/2(l + M 2 )] 

(C7) 



where a « 1 — 2 is a parameter approximating the complex 
shape of the collapse barrier. 

The left-hand panel in Figure 



CI 



compares 



and 



/con as a function of Q(disc) for several Mach numbers. 
For Q(disc) < 1, the predicted collapsed fraction is ~ 1, 
as expected since the disc is unstable in an average sense. 
For Q(disc) > 1, GMCs only form where turbulent den- 
sity fluctuations are large enough to locally bring Q(h) < 1 



15 Since h = [ct / (V^4>o')]i" (eq. fTol , this is equivalent to 
parametrizing the disc by a constant c^(h) and viewing gravi- 
tationally collapsed regions as corresponding to upward density 
fluctuations about the mean p — cr relation implied by equation 
| |C2| |. In reality, both CT(h) and p(h) are likely to fluctuate spa- 
tially, possibly with non-trivial correlations. These effects are not 
captured by our analytic estimate. 



(in the Toomre heuristic picture). Since density fluctuations 
are suppressed by rotation and finite mass effects on scale 
~ h, the collapse fraction decreases rapidly as Q(disc) is in- 
creased above unity. The right-hand panel shows power-law 
approximations used for analytic estimates in |3.3| values 
of the exponent a ~ 3 — 5 (see eq. |31[ ) provide reasonable 
approximations for Mach numbers M = 10 — 100. 
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